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ABSTRACT 

The so-called star-forming main sequence of galaxies is the apparent tight relationship 
between the star formation rate and stellar mass of a galaxy. Many studies exclude 
galaxies which are not strictly ‘star forming’ from the main sequence, because they do 
not lie on the same tight relation. Using local galaxies in the Sloan Digital Sky Survey 
we have classified galaxies according to their emission line ratios, and studied their 
location on the star formation rate - stellar mass plane. We find that galaxies form 
a sequence from the ‘blue cloud’ galaxies which are actively forming stars, through 
a combination of composite, Seyfert, and low-ionization nuclear emission-line region 
galaxies, ending as ‘red-and-dead’ galaxies. The sequence supports an evolutionary 
pathway for galaxies in which star formation quenching by active galactic nuclei (AGN) 
plays a key role. 
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1 INTRODUCTION 

Galaxies can be separated into two distinct classes: those 
that are actively forming stars and are blue in colour, and 
those that have no appreciable star formation and are red in 
colour. Current consensus is that galaxies start their lives in 
the active ‘ blue’ category, and then evolve to be “red-and - 
dead” (e.g. iGoncalves et al.l l2012l . iMoustakas et al.l l2013l b 
Evolution between these two classes must involve processes 
which quench the star formation, and quench it relatively 
quickly (on the order of 1 Gyr), because galaxies with in¬ 
termediate propertie^loca^d in the ‘green valley’, are not 
often seen (|Faber et al.l . 1200^ 1. The dearth of observed green 
valley galaxies could also be explained if they become ob¬ 
scured, or if their number density is underestimated due to 
the sho rt variability/duty c ycle time-scales of black hole ac¬ 
cretion (jHickox et al.L 12014 ). A popular candidate quenching 
mechanism is negative feedback from active galactic nuclei 
(AGN), because AGN host galaxies tend t o have green or 
red optical colours (e.g. iNandra et al]l2007l b 

Blue star-f orming galaxies form a tight (~ 0.2 dex in¬ 
trinsic scatter; ISDeagle~ 'Zl l20i4) sequence between star 
formation rate (SFR) and stellar mass (M*) often re- 
ferred to as the galaxy ‘main sequence’ of star formatio n 
(jNoeske et all . l2007l . IPaddi et all . [20071 . lElbaz et~M] . l2007l ) . 
The main sequence (MS) is generally characterized by a sin- 
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gle power law, often expressed as 

log(SER) =/31og(M.) + a (1) 

with t ypical values for the slope jd — 0.7— l.O. iSpeagle et al.l 
(|2014l ) found that the slope, /3, and normalization, or inter¬ 
cept, a, increase and decrease respectively, as a function of 
the age of the Universe t (in Gyr): 

/3 = (0.84 ±0.02 - 0.026 ± 0.003 xt), (2) 

a = (-6.51 ±0.24 ±0.11 ±0.03 X t). (3) 

SER at 10^° M© evolves from ~10 to ~1 M© yr“^ from 
z~l to the present day. The existence of the MS implies 
that star formation in the majority of galaxies is governed 
by quasi-steady processes, and may imply that only a small 
fraction of local galaxies undergo more c haotic processes 
such as major mergers (|Lotz et al.l . l201lh . Galaxies with 
quasi-exponentially increasing SERs and stellar masses at 
high redshif t will form an MS relation with M* oc SER 
(|Davd . lioOSi ). The fitted parameters of the MS relation re¬ 
veal even more about the nature of star formation and stel¬ 
lar mass build up: the slope of the MS is related to the 
low mass slope of the mass function of star-forming galax¬ 
ies (jPeng et al.l . 120101 ). and the evolution of the MS nor¬ 
malization is tho ught to result from evolution in gas densi¬ 
ties with redshift (IWhitaker et al.l . l2012l . iMagdis et al.l . l2012l . 
iSoeagle et al.l . l2014l ). 

Different studies report incon sistent slopes and shape s 
of the MS. F or exam p le, studies bviProrv Alvar"^ (|2008|). 
[Karim et al ] IWhitaker et alJ ll2ni4h fc II ee et al J 
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(|2015l l suggest that the MS requires a more com plex, stel¬ 
lar ru ass dependent fit. iKarim et al ] (I 2011 I 1 & Ihee et ahl 
(|2015l l report a flattening of the MS at stellar masses 
M* > On the other hand, after including both 

star-forming and quiesc ent galaxies from the Sloan Digi¬ 
tal Sk y Survey (SDSS; I York et ^ l2000l i . iRenzini Pend 
(|2015l i found no deviation from a simple power-law MS re¬ 
lation in local galaxies, which suggests that the deviations 
might originate from inconsistent definitions of star-forming 
galaxies and/or selecti on effects. However, the analysis of 
IRenzini fc Pend (|2Q15l i excluded galaxies exhibiting AGN 
activity, because it has been historically difficult to reliably 
measure their underlying SFRs. 

Studies which include AGN in colour-mass diagrams 
have shown that AGN host galaxies pop ulate the red- 
seque nc e or the transit io nal gr ee n valley ([Nandr a et , 
2 OO 7 I. ISchawinski et al.L l2007l . ISmolcic et al.L 12009 . 
Alatalo et al.Ll2014bl b The inclusion of AGN in colour-mass 
diagrams thus allowed the role AGN feedback in the 
quenching th e star formation of ea rly type galaxies to be 
inferred fe.g. ISchawinski et ahllioo^ . 

In this Letter we use the SDSS to investigate the effects 
of including all spectral classes of galaxies, in the main se¬ 
quence analysis. We will show that AGN host galaxies also 
populate intermediate regions on the M* - SFR plane. 


2 DATA SAMPLE 

The dat a used in this study com es from the SDSS data re¬ 
lease 7 (|Abazaiian et al.l . l2009l i . with all properties calcu¬ 
lated by the MPA/JHU groupj. We only use data for galax¬ 
ies with redshifts greater than z > 0.04 to ensure that the 
SDSS fibre covers at least 30% (median cover age of 38%) of 
the t ypical galaxy to avoid aperture effects (jKewlev et al.l . 
l2005l b We also constrain our sample to z<0.1 t o avoid incom- 

S iess and significant evolutionary effects (jKewlev et al.L 
. This redshift restriction results in a sample of 273,274 
galaxies with MPA/JHU derived emission line fluxes, SFRs 
and stellar masses. We further restrict our sample by ap¬ 
plying a signal to noise limit of 3 on the Ha emission line, 
resulting in a sample of 202,708 galaxies. We have also re¬ 
trieved reddening corrected rest-frame photometry for these 
galaxies from the SDSS DR 7 data base. We use the mag¬ 
nitudes resulting from fitting either a de Vaucouleurs pro¬ 
file or a pure exponential profile, which have been cor- 
re cted for galact i c exti nction following the determinations 
of ISchlegel et al ] (Il998ll . 

We use the SFRs from the MPA/JH U catalogue, which 
were c alculated based on the methods of iBrinchmann et al.l 
(|2004| . with modifications to the methods for non-SF galax¬ 
ies. For galaxies classified as SF, the SFR is cal culated by 
the M PA/JHU team using the Ha calibration of iKennicuttI 
They use th e Balmer decrement and method of 
Ig larlot Falll (l200Ql i to correct for extinction and assume 
a Kroupal (|200lh initial mass function (IMF). For galaxies 
falling outside of the SF category, the MPA group calculated 
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^ Note that using the Kroupa IMF, iKennicutt fc EvaTiS (l2012li 
present a Ha SF R which is 0.17 dex lower than the calibration in 
iKennicuttI lll998l) . 


their SFRs estim ated based on the stre ngth of the D4QQQ 
break (defined in iBalogh et al.lll998l L In iBrinchmann et al.l 
||20 o 3) , an empirical relationship was derived between Ha 
SFR and D4000, for SF galaxies. An updated version of the 
resulting relationship is used to calculate SFRs for galaxies 
in which the Ha flux may not be an accurate indicator of 
SFR (e.g. in the presence of AGNs, or quiescent galaxies 
with l ow S/N emission lines). Unlike in IBrinchmann et al.l 
(|2004I L the aperture corrections were made for SFRs by the 
MPA/JHU team by f itting the photome try of the outer re¬ 
gions of the galaxies (| Salim et al.l . l2007l i . 

Stellar masses are calculated by the MPA/JHU team us¬ 
ing a fit to photometry. The stellar masses in the MPA/ JHU 
catalogue, h ave been found to be consistent with other es¬ 
timates (see iTavlor et al.ll201ll . Ichang et al.ll2015l b In the 
following analysis we use the median SFR and M* from the 
probability distributions for each galaxy, however using the 
mean does not change our results. 


2.1 Classification 

We c lassify each galaxy according to the iKewlev et al.l 
(|2006l i scheme which makes use of the emission line ra¬ 
tios [NH]/Ha, [SH]/Ha, [OI]/Ha, and [OHI]/H/3. These 
em ission line rati o s wer e pro posed in diagnostic di a grams 
by iBaldwin et al.1 (|l98lh and fVeilleux OsterbrockI (|l987l i 
to classify the dominant energy source in an emission-line 
galaxy. We refer to the optical classification diagrams as the 
‘standard optical diagnostic diagrams’. 

We find that out of the 202,708 galaxies in our sample, 
60.4% are classified as SF, 12.2% are classified as compos¬ 
ite, 4.1% are classified as Seyferts , 6.5% are classified as 
LINERs (low-ionization nuclear emission-line regions) and 
17.0% are classified as ambiguous- i.e. they do not fall into 
the same category on all the three diagnostic diagrams. The 
fraction of ambiguous galaxies becomes smaller if signal-to- 
noise limits are applied to the forbidden lines. Most galaxies 
classified as ambiguous lie close to the Seyfert 2-LINER or 
composite-AGN boundaries. We expect some of these am¬ 
biguous galaxies to have the properties of composite, Seyfert 
2, or LINER galaxies. 


3 THE ‘GALAXY MAIN SEQUENCE’ 
INCLUDING AGN HOSTS 

The canonical main sequence is shown for star-forming 
galaxies in the f irst panel of Eigur eU The local SDSS rela¬ 
tion derived bv lElbaz et al.l (|2007l i from galaxies with blue 
optical colours is represented by the black dashed line in all 
panels of Eigure[T] The SE galaxies form a relationship be¬ 
tween log(MH<) and log(SER), with most galaxie s class ified 
as SE lying along the MS relation of lElbaz et al.l (|2007l L 
Unlike the pure SE MS, there is almost no relation be¬ 
tween stellar mass and SER for composite galaxies, other 
than the fact that composite galaxies all have high stellar 
masses This stellar mass limit may also exist 

for quasar hosts (jMatsuoka et al.L l2014l i . Gomposite line ra¬ 
tios may be ca used by a combin ation of star formation and 
AGN activity (jYuan et al.Ll201C)h. a comb i nation of star for¬ 
matio n and shock excitation ( Allen et a D. l2008l.lRich 


or a combination of all three processes ([Leslie et al 
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Table 1. For each class: the percentage of the total sample 
(202708 galaxies), median signal to noise in the Ha emission, 
median (log) stellar masses and star formation rates and their 
corresponding inter-quartile range. 


Class 

% 

S/N 

log(M*) 

\og{SFR) 

SF 

60.4 

24 

10.02 

0.65 

0.09 

0.51 

Composite 

14 

35 

10.56 

0.44 

0.04 

0.70 

Seyfert 2 

4.1 

15 

10.58 

0.51 

-0.45 

0.93 

LINER 

6.5 

5 

10.74 

0.46 

-0.79 

0.53 

Ambiguous 

17.0 

6 

10.65 

0.64 

-0.53 

0.79 


I2OI4I ). The contours of composite galaxies extend above the 
top star-forming MS. This region of high SFR could cor¬ 
respond to gala xies which are unde r going intense bursts of 
star formation (|Schiminovich et~al] . 12007 ). The composite 
sequence extends from the massive end of the main sequence 
down to the region where red galaxies, with little active star 
formation lie. The mean stellar mass of composite galaxies is 
higher than that of pure SF galaxies, however the mean SFR 
is lower. This could imply that suppression of star formation 
has commenced by the composite stage. 

Further bridging the gap between star-forming galaxies 
and quiescent galaxies are the Seyfert 2s. The majority of 
Seyfert 2s lie in the quiescent region of the diagram. 

LINERs do not fall on the star-forming MS, but lie 
below it. Most LINERs in our sample clearly contain very 
little current star formation activity. On average, LINERs 
have a higher stellar mass than Seyfert galaxies (median 
log(M*) of 10.74 an d 10.58 Me:) r e specti vely), consistent 
with the findings of iKewlev et al.l (|2006l L The strongest 
difference between Seyferts and LINERs is their Edding¬ 
ton ratios. Nuclear LINERs that contain A GN are found 
to ac c rete less efficiently t han Seyfert AGN (jKewlev et al.l . 
I 2 OO 6 I ) . iKewlev et al.l (| 20061 ) also found that LINERs are gen¬ 
erally found in galaxies with an older stellar population than 
Seyfert hosts. 


LIN ER-like emis s ion may also be caused by 
shocks (|Kewlev et al.L l200lh . Shocks can be pr o duced 
by AGN- r e lated processes (e.g. IPavis et~M ] 
lAlatalo et all 2014alla) or in star b urst driven winds 
(e.jShar;) Bland-HawthornI I 2 OIOI . I Rich et al. I HH, 
IHo et, aLl2ni4l l. Winds driven by a powerful starburst, may 
also drive gas out of a galaxy disc and thus suppress the 
SER. 


The final panel of Eigure[T] shows all the galaxy classes 
together, with contours at 10, 30, and 60% of the maximum 
number density for each class of galaxies. In general compos¬ 
ite galaxies have a higher SER and lower stellar mass than 
Seyfert2s, which in turn have a higher SER and lower stellar 
mass than LINERs, with ambiguous galaxies having the low¬ 
est SERs (see Table 1). As we will discuss in Section 4, this 
is consistent with an evolutionary sequence through which 
galaxies evolve from main sequence star-forming galaxies 
to being quiescent galaxies through the composite, Seyfert 
and/or LINER classes. 

The colour-mass diagram in Eigure[2l showing u-r as a 
function of log(M*) reinforces the story of the MS in Eig- 
ure [U Purely star-forming galaxies have blue colours and 
span a range of stellar masses. Gomposite galaxies have red¬ 
der colours and are more massive than purely star-forming 
galaxies. On average Seyfert 2s have redder colours than 


composite galaxies and are more massive, however there 
is a significant arnount of overlap between the two classes. 
ISchawinski et al.] (|2007l ) found similar results for a sample 
of early type galaxies and suggested that the overlap im¬ 
plies that galaxies may transition between t he Seyfert and 
compo site regions. fSchawinski et al.l (|2007l ) fc ISmolcic et al.l 
(I2OO9I ) report a similar trend in the colour-mass diagram, 
with a sequence existing between the blue cloud and the 
red sequence; from blue star-forming galaxies at the lowest 
mass, through Seyferts in the green valley, to LINERs at the 
highest masses and reddest colours. 


4 EVOLUTIONARY PATHWAYS 


We interpret the trends seen in Eigures [T] and [2] as evi¬ 
dence that optical emission line classes provide important 
pathways for galaxies reaching the red sequence. Galax¬ 
ies start their lives on the main sequence forming stars 
and gaining stellar mass at a steady rate until they reach 
M* ~ The transition away from the main sequence 

for our sample only occurs at stellar masses > lO^^M©. 
At these high stellar masses mass-quenching is expected 
to dominate over environmental quenching, which could in- 
yolve strangulatio n, ram-pressure stripping or harassment 
(|Peng et al.l I 2 OIOI ). Our work indicates that galaxies may 
transition from star-forming through one of, or a combi¬ 
nation of composite, Seyfert, and/ o r LIN E R stages before 
achiev ing quiescenc e. Karim et al. I (|201lh . I Whitaker et al.l 
(|2014l ) & iLee et al.] (|2015l ) report a turnover in the MS at 
stellar masses ~ 10^°. The inclusion of composite and AGN 
galaxies in these high redshift studies may be the major fac¬ 
tor causing the flattening of the main sequence at stellar 
masses > lO^^M©. 

The loci of SER vs M* for the composite, Seyfert, and 
LINER classes are almost perpendicular to the star forma¬ 
tion main sequence showing decreasing sSER (specific star 
formation rate = S ER/M^) for ea c h clas s consistent with 
recent results f rom Sh imizu et~al] (|2015l ) and with semi¬ 
nal results from ISalim et ahh 20071 ) . The SERs of composite 
galaxies lie above and below the MS relation. The composite 
galaxies a bove the MS are likely sta rburst galaxies similar 
to those in ISchiminovich et ^ (|2007l ) . Quenching may have 
commenced in composite galaxies below the MS. Gomposite 
galaxy quenching could be tied to large scale gas outflows 
from winds driven by bursts of star formation, with the com¬ 
posite spectra resulting from shocks embedded in the out¬ 
flows. Alternatively, the quenching could be tied to AGN ac¬ 
tivity with the composite spectra resulting from a mixture 
of star formation and AGN emission. AGN jets can drive a 
wind which is able to blow out gas from a galaxy. Activity 
from an AGN can also heat gas, preventing star formation. 
Galaxies may evolve off the MS towards the red-sequence, 
via one or more of the composite, Seyfert and LINER stages. 

We emphasize that different phenomena may produce 
the same emission line ratios for galaxies within each class. 
Gomposite emission spectra are not only produced by SE 
and AGN, but are often shock excited. LINER-like emis¬ 
sion may be caused b y pos t AGB stars (|Singh et al.L l2013l . 
iGid Eernandes et ah. I 2 OO 9 I) or shocks, a s well as by low- 
luminosity AGN (|Ho Eilippenkj . Il993l ) . To elucidate the 
nature of galaxies in these classes, and their relationship to 
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Figure 1. SFR as a function of stellar mass for each class of galaxies. Stellar mass and SFRs are in units of Mq and MQyr“^ respectively 
throughout this work. The top-left to bottom-right panels include galaxies classified as purely star-forming, composite, Seyfert 2, LINER, 
or Ambiguous, and the final panel contains all classe s together. A black dashed line in each panel represents the local MS relation for 
blue SDSS galaxies determined bv lElbaz et al.l (l2007l b Contours and colours represent the number density of galaxies in a single class 
only. Dividing the SFR-M* space into 150x150 bins, contours are drawn at 10, 30, and 60% of the maximum number density. The blue 
contours of the star-forming galaxies are included in all panels to indicate the location of the star-forming MS. Contours from previous 
panels are also shown in the final panel, which displays the MS of all galaxies. 



8.0 8.5 9.0 9.5 10.0 10.5 11.0 11.5 12.0 8.0 8.5 9.0 9.5 10.0 10.5 11.0 11.5 12.0 8.0 8.5 9.0 9.5 10.0 10.5 11.0 11.5 12.0 



Figure 2. Galaxy reddening-corrected rest-frame colour as a function of stellar mass for each class. The red contour in the Pure SF 
panel corresponds to the same galaxies in the red contour of Figure 1. Contours are drawn at 10, 30 and 60% of the maximum density 
of each class of galaxy. 


star formation quenching, our future work includes wide in- 5 CONCLUSION 
tegral field spectroscopy of statistically significant datasets 
in each class. 

Using star formation rates and stellar masses from the 
MPA/JHU measurements of SDSS galaxies, we have in¬ 
cluded composite and AGN host galaxies on the local star- 
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forming main sequence and colour-mass diagrams. Our re¬ 
sults suggest that: 

• AGN activity (composites, Seyferts, LINERs) appears 
to play a very important role in “quenching” star formation 
in massive (M* > 10^° M©) galaxies; 

• the loci of SFR vs. M* for composites, Seyferts and 
LINERS are nearly perpendicular to the SE MS. Compos¬ 
ites show the broadest range in SER (+0.8 to -1.8 dex 
above/below the MS), while most Seyferts lie below the MS 
{SFR ~ 0 to -2 dex), and LINERs lower still ( SFR ~ -1 
to -2 dex); 

• composites, Seyferts and LINERs plausibly form a set 
of evolutionary pathways by which massive SE MS galaxies 
transition into a “red and dead” quiescent phase. 

Euture work with large spatially resolved inte gral field unit 
data sets, such as those provided by M aNGA ([Bundv et al.l . 
I 2 OI 5 I ) and SAMI (jCroom et al.l . I 2 OI 2 I I , will likely shed more 
light on this sequence, by allowing the study of the spatial 
distribution of star formation in galaxies as a function of 
galaxy properties such as morphology, stellar mass, BPT 
classification, and environment. 
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